Proton inelastic diffraction by a black nucleus and the size of excited nuclei 
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We systematically derive a length scale characterizing the size of a low-lying, /3 stable nucleus 
from empirical data for the diffraction peak angle in the proton inelastic differential cross section of 
incident energy of ~ 1 GeV. In doing so, we assume that the target nucleus in the ground state is 
a completely absorptive "black" sphere of radius a. The cross section iva 2 , where a is determined 
in such a way as to reproduce the empirical proton diffraction peak angle in the elastic channel, is 
known to agree with empirical total reaction cross sections for incident protons to within error bars. 
By comparing the inelastic diffraction patterns obtained in the Fraunhofer approximation with the 
experimental ones, one can likewise derive the black sphere radius ai for the excited state with spin 
I. We find that for 12 C, 58 > 60 > 62 . 64 Ni, and 208 Pb, the value of ai obtained from the inelastic channel is 
generally larger than the value of a from the elastic channel and tends to increase with the excitation 
energy. This increase is remarkable for the Hoyle state. Finally, we discuss the relation between ai 
and the size of excited nuclei. 

PACS numbers: 21.10.Gv, 24.10.Ht, 25.40.Ep 



I. INTRODUCTION 

The size of atomic nuclei is one of the most fundamental quantities that characterize matter in the nuclei. It is 
well known for /3 stable nuclei in the ground state thanks to systematic measurements of electron and proton elastic 
differential cross sections [l| . This helps clarify the equation of state of nuclear matter near the saturation point 0] . 
For excited states of (3 stable nuclei, however, it is not straightforward to deduce the nuclear size, because elastic 
scattering off an excited target is hard to measure. Alternatively, one can use proton inelastic differential cross sections 
in deducing the size of excited nuclei, but all one may know is the transition density, which only implicitly reflects 
the density distribution of the excited nuclei. 

Although one can in principle determine the charge and matter density distributions of a target nucleus from 
measured electron and proton elastic differential cross sections, it is practically time-consuming even if one uses 
optical potential models and Glauber's multiple scattering theory in an approximate manner. Instead of sticking to 
such microscopic derivations of the nuclear size, in Ref. [3[ we constructed a phenomenological method for deducing 
the nuclear size by focusing on the peak angle in the proton-nucleus clastic differential cross section measured at 
proton incident energy T p ^800-1000 MeV, where the corresponding optical potential is strongly absorptive. In this 
method, we regard a nucleus as a "black" (i.e., purely absorptive) sphere of radius a, and we determine a in such a 
way that the first peak angle of the Fraunhofer diffraction by a circular black disk of radius a agrees with that of 
the measured diffraction. This method is reasonable as long as the scattering is close to the limit of the geometrical 
optics. This condition is fairly well satisfied at least for T p > 800 MeV, where the wave length of the incident proton 
is sufficiently shorter than a even for 4 He. The black sphere picture is originally expected to give a decent description 
of total reaction cross sections for any kind of incident particle that tends to be attenuated in the nuclear interiors. In 
fact we showed that for proton beams incident on stable nuclei, the cross section of the black sphere of radius a thus 
determined is consistent with the measured total reaction cross section j4j. If we multiply a by y3/5, furthermore, 
a ratio between the root-mean-square and squared off radii for a rectangular distribution, the result for stable nuclei 
of A > 50 shows an excellent agreement with the root- mean-square radius, r m , of the matter density distribution 
as determined from conventional scattering theories so as to reproduce the overall diffraction pattern and analyzing 
power in the proton elastic scattering 

In this paper, we apply the black sphere picture to analyses of proton inelastic scattering data for T p ^1000 MeV. 
Basically, this application is a straightforward extension of the case of proton elastic scattering, which is closely related 
to the method developed by Blair |5j for alpha scattering by assuming elastic diffraction by a circular black disk of 
radius Rq and inelastic diffraction by a black nucleus with small multipolar deformations from a sphere of radius Rq. 
The present extension is, however, accompanied by a nontrivial choice of the inelastic diffraction peak whose angle is 
to be fitted to the empirical value and by a variation of the black sphere radius from the value determined from the 
elastic diffraction peak angle. The resulting black sphere radius does not correspond to the size of the nucleus excited 
by the incident proton, but rather is related to the transition density and thus expected to lie between the sizes in the 
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ground state and in the excited state. Even so, as we shall see, systematic derivation of the black sphere radius from 
the inelastic channels is useful for predicting how the size in the excited state depends on the excitation energy E ex . 

In Sec. II we extend the black sphere approach developed for proton elastic diffraction to the case of proton inelastic 
diffraction. The results for the black sphere radii are illustrated in Sec. III. 



II. BLACK SPHERE APPROACH 

We begin by summarizing the black sphere approach to proton elastic diffraction Q. The center-of-mass (cm.) 
scattering angle for proton elastic scattering is generally given by c . m . = 2sin~ 1 (q/2p) with the momentum transfer 
q and the proton incident momentum in the cm. frame p. For the proton diffraction by a circular black disk of radius 
a, we can calculate the value of 6 c . m . at the first peak as a function of a. (Here we define the zeroth peak as that 
whose angle corresponds to c .m. = 0.) We determine a in such a way that this value of # c . m . agrees with the first 
peak angle for the measured diffraction in proton- nucleus elastic scattering, 9m- The radius, a, and the angle, 8m, 
are then related by 



2pasin(0 M /2) = 5.1356- 



(1) 



By setting 



r BS = 7375a, (2) 

we found Q that at T p > 800 MeV, res, estimated for heavy stable nuclei of A > 50, is within error bars consistent 
with the root- mean-square nuclear matter radius, r rn , deduced from elaborate analyses based on conventional scat- 
tering theory. Thus, expression © works as a "radius formula." The factor -\/3/5 comes from the assumption that 
the nucleon distribution is rectangular; the root-mean-square radius of a rectangular distribution is a cutoff radius 
multiplied by \/3/5. For stable nuclei with A < 50, however, the values of res are systematically smaller than those 
of r m 4]. The scale a is nevertheless meaningful because the values of na 2 for C, Sn, and Pb agree well with the 
proton- nucleus reaction cross section data for T p > 800 MeV [H . This indicates that a can be regarded as a "reaction 
radius," inside which the reaction with incident protons occurs. In a real nucleus, this radius corresponds to the 
radius at which the mean free path of incident protons is of the order of the length of the penetration. We remark 
that even for deformed nuclei, this interpretation works well unless the degree of deformations is extremely large. 

We now proceed to generalize the black sphere picture to the case of proton inelastic scattering by following a 
line of argument of Blair We assume that the final low- lying spin-Z excitation of a target even-even nucleus is 
characterized by small multipolar deformations of a black sphere of radius a;. As we will see, this radius generally 
differs from the value of a obtained from the measured peak angle of elastic diffraction off the same target, as well as 
among different levels with the same I. Therefore, the radius a; is conceptually different from the radius Rq used in 
Ref. [5]. In the adiabatic approximation, which is valid for the proton incident energies. T p > 800 MeV. of interest 
here, we may set the reaction radius as 



Rl(6,<p)=ai 
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ai m Yi m (9, 4>) 



(3) 



where 9 is measured with respect to the incident beam axis, and ai m are the deformation parameters. For given ai m 
and q, we write the scattering amplitude in the cm. frame of the proton of initial momentum p and the black nucleus 



f(q;a hn ) = — 
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(4) 



where b is the impact parameter perpendicular to p. Here we assume that the final proton momentum p + q has the 
magnitude equal to that of the initial momentum p because T p is far larger than the excitation energy _B ex . Then, 
the cm. scattering angle again becomes 9 c . m , — 2 sin~ 1 (g/2p). With the convention that <fi is measured with respect 
to the projected final proton momentum onto a plane perpendicular to p, the scattering amplitude up to linear order 
in ai m reads 



/(q; ai m ) 



ip 
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-iqb cos 4> Mh 



(•5) 
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The first integral would give the amplitude for elastic scattering off the excited nucleus with spin I: 

ipafJi(2pai sin(fl c , m ,/2)) 
2pai sin(0 c .m./2) 



(6) 



Note that this amplitude is independent of ai rn . Anyway, the first integral is beyond the scope of real experiments. 
The second integral in Eq. ([5]) leads to the amplitude for inelastic scattering that excites a target nucleus in the 
ground state to the final state of spin I 
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x on m J\ m \{2pai sin(0 c . m ./2)), (7) 
where (2ro)!! = 2 • 4 • 6 • • • 2n. The corresponding differential cross section is 

da ' {l-m)\{l + m)\ 

dfl^ > 00 2^ \(l - m )\\(l + m )W 

m=-l,-l+2,-- 

xaL4 l |(2pa/sin(0 c . n ,/2)), (8) 

which gives the inelastic Fraunhofer diffraction pattern. 

We next determine ai by comparing the peak angle of the black sphere inelastic diffraction as described by Eq. 
([8|) with the corresponding empirical value 6 mi ■ Generally, the inelastic diffraction pattern at forward angles is more 
complicated than the elastic one. We thus avoid the peak that is the nearest to 9 C m — and consider which peak to 
be chosen among the rest. The guiding principle for this choice is simply to search for the peak that systematically 
corresponds to the first peak in the elastic diffraction pattern just like the case of electron diffraction [6J . We thus 
obtain 

2pa o sm(0Mo/2) = 7.015 (9) 
2pa 2 sm(0M2/2) = 6.783- •-, (10) 
2pa 3 sm(6 M3 /2) = 8.209- ••, (11) 

2pa 4 sin(6> M4 /2) = 9.617 • • • . (12) 

The values of the right side in Eqs. (|9|)- ([T2j) correspond to one of the values of 2pai sin(6> c . m ./2) = x where Eq. (|SJ) 
is locally maximal. In fact, the values of x are as follows: For I = 0, x = 0, 3.831- • -, 7.015- for I = 2, x = 0, 

3.251- • -, 6.783- for 1 = 3, x = 1.982 • ■ -, 4.605- • -, 8.209- • -, . . .; for I = 4, x = 0, 3.675- • -, 5.852- • -, 9.617- • -, 

. . .. The choice of other diffraction peaks in extracting 9 mi would in many cases produce a; that is significantly 
different from a. We remark that the avoidance of the local diffraction maxima at forward angles justifies the neglect 
of Coulomb effects in the present black sphere approach [Bj]. 

The extraction of 9 mi from the measured inelastic differential cross sections is not always straightforward for several 
reasons. Firstly, in the case in which the empirical data lack such diffraction maxima as are supposed to occur in 
the Fraunhofer diffraction, counting the empirical diffraction maxima does not work. It is thus indispensable to 
compare the overall measured diffraction pattern with the elastic one whose first peak is easily distinguishable, before 
identifying the local maxima that give 9 mi- Secondly, even after the local maxima are successfully identified, we need 
to allow for various uncertainties accompanying determination of 9 mi ■ 

As in the case of elastic scattering |3[, we basically determine the values of 9 mi from the scattering angle that 
gives the maximum value of the cross section among discrete data near the identified diffraction maximum. Some 
examples are shown in Fig. 1. When the data stagger in such a way that the scattering angle that gives the local 
maximum is significantly away from the diffraction peak position deduced from the overall plot of the data, we select 
the data point that obviously seems the closest to the peak position. Such determination of 9mi is accompanied by 
uncertainties in the measured scattering angle and systematic errors that are dependent on the way of deducing the 
peak position. The uncertainties in the measured angle, which are due mainly to the absolute angle calibration, are 
typically of order or smaller than ±0.03 deg Q for existing data for proton scattering off stable nuclei. On the other 
hand, the systematic errors can be estimated by assuming that the true peak is located in the region enclosed by 
the two data points that are the closest neighbors of the selected data point. The systematic errors thus estimated 
dominate the error bars of 9 mi, and hence we will ignore the uncertainties in the measured angle. 
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FIG. 1: (Color online) The inelastic differential cross section calculated from the Fraunhofer diffraction formula (|8]) for p- 12 C 
and p- 208 Pb (T p = 1040 MeV) as well as p- 58 Ni (T p = 1047 MeV). The experimental data (large dots) are taken from Refs. 
[1)01; f° r comparison, the elastic data (small dots) are also plotted. In each panel, the vertical lines denote the peak angle 8 mi 
(including error bars) at which we fit the calculated peak angle to the empirical one, and the normalization of Eq. (|8]) is set in 
such a way that the value of the differential cross section at 8mi agrees with the empirical one. 



III. BLACK SPHERE RADII 



We finally obtain ai from the determined 9 mi via Eqs. ([9l)-(fl~2]). For comparison, we likewise determine Om and 
then a via Eq. ([!} . We remark that a part of the errors of ai and a that come from uncertainties in the proton incident 
energy, which are typically a few MeV , are negligible. The results for a and aLpbtained from empirical scattering 
data off 12 C, 58,60,62,64 Nij and 208p b &t proton incident energy of about 1 GeV Hill are listed in Table I. In collecting 
the data, we have made access to the Experimental Nuclear Reaction Data File (EXFOR) [l(| ■ In the absence of the 
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TABLE I: Black sphere radii obtained from proton elastic and inelastic scattering data. 



Target 


Final state 


S cx (McV) T P 


(MeV) 


a or ai (fm) 


12 C 


g.s. 





1040 


2.75±0.06 




2+ 


4.44 


1040 


2.70±0.06 




0+ 


7.65 


1040 


3.20±0.07 


58 Ni 


g.s. 





1047 


4.79±0.18 




2+ 


1.45 


1047 


4.91±0.14 




4+ 


2.46 


1047 


5.22±0.11 




3" 


4.47 


1047 


5.09±0.13 


60 Ni 


g.s. 





1047 


4.78±0.18 




2+ 


1.33 


1047 


4.91±0.14 




4+ 


2.50 


1047 


5.22±0.12 




3" 


4.04 


1047 


5.09±0.13 


62 Ni 


g.s. 





1047 


4.96±0.19 




2+ 


1.17 


1047 


5.05±0.15 




4+ 


2.34 


1047 


5.22±0.12 




3~ 


3.76 


1047 


5 22+0 1 3 

'.1.1—1— 1 V / . 1 « J 


64 Ni 


g.s. 





1047 


4.96±0.19 




2+ 


1.35 


1047 


5.04±0.15 




4+ 


2.61 


1047 


5.22±0.11 




3" 


3.55 


1047 


5.09±0.13 


208p b 


g.s. 





1040 


7.49±0.35 




3" 


2.62 


1040 


7.29±0.35 



inelastic scattering data for the 4 + state of 12 C and the 2 + and 4 + states of 208 Pb, the corresponding black sphere 
radii are unavailable. Note that 03 is unavailable for 12 C despite the presence of the data for the 3~ state of 12 C, 
because the corresponding peak is missing in the measured differential cross section. 

To see the systematic behavior of the black sphere radii, we plot the radii as function of E ex in Fig. 2. We find 
that the black sphere radius tends to increase with E cx , with a few exceptions in which case the black sphere radius 
decreases with E cx in its central value but can be regarded as unchanged allowing for the error bars. This is consistent 
with the behavior of the transition radii obtained systematically from electron inelastic scattering off 208 Pb Q. 

Recall that the black sphere radius a; is related to the transition density and thus expected to lie between the 
sizes in the ground state and in the excited state. It is thus reasonable that the differences between a and ai for the 
low-lying excited states of interest here are generally small. The only exception that we discovered is the 12 C + state 
(the Hoyle state), for which a^/a ~ 1.16. The nucleus in the Hoyle state is thus expected to be larger than that in 
the ground state by more than 16 %, a feature consistent with the a-clustering picture of the Hoyle state fi"Tl |. 

In summary, we generalized the black-sphere method for deducing the size of the ground-state nuclei from proton- 
nucleus elastic scattering data to the case of proton-nucleus inelastic scattering data and provided implications for 
the size of excited nuclei. In the present analysis, we confined ourselves to even-even nuclei. Extension to odd nuclei 
would be straightforward if the conventional collective model works [f| . It is also important to note that the validity of 
the inelastic Fraunhofer diffraction formula used here tends to decrease with increasing scattering angle. We assume 
that the peak selected for determination of ai, as shown in Fig. 1, is in the valid regime, which could be checked by 
microscopically clarifying a relation between the transition density and ai (l2l |. We hope that the present analysis 
could develop into a systematic drawing of the black-sphere radii of isomers and nuclei in other characteristic excited 
states over a chart of the nuclides. 
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FIG. 2: (Color online) The black sphere radius as a function of the excitation energy. The target nucleus is 12 C, 58 > 6 °. 62 > 64 Ni, 
and 208 Pb. For Ni isotopes, the squares (with line), circles (with short-dashed line), triangles (with long-dashed line), and 
crosses (with dash-dotted line) denote the results for 58 Ni, 60 Ni, 62 Ni, and 64 Ni, respectively. 
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